Multiple quantum phases in artificial double-dot molecules 
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We study coupled semiconductor quantum dots theoreti- 
cally through a generalized Hubbard approach, where intra- 
and inter-dot Coulomb correlation, as well as tunneling ef- 
fects are described on the basis of realistic electron wave- 
functions. We find that the ground-state configuration of 
vertically-coupled double dots undergoes non-trivial quantum 
transitions as a function of the inter-dot distance d; at inter- 
mediate values of d we predict a new phase that should be 
observable in the addition spectra and in the magnetization 
changes. 



Semiconductor quantmii dots (QDs) are nano- or 
mesoscopic structures that can be regarded as 'artificial 
atoms' because of the three-dimensional carrier confine- 
ment and the resulting discrete energy spectrum. They 
are currently receiving great attention because they can 
be designed to study and exploit new physical phenom- 
ena: On the one hand, the nature and scale of electronic 
confinement allows the exploration of regimes that are 
not accessible in conventional atomic physics; on the 
other hand, they lead to novel devices dominated by 
single- or few-electron effects [|l|,|[ . 

One of the important challenges at this point is to un- 
derstand the fundamental properties of coupled quantum 
dots, the simplest structures that display the interactions 
controlling potential quantum-computing devices. Here 
also the interdot coupling can be tuned through exter- 
nal parameters, far out of the regimes known in 'natural 
molecules' where the ground-state interatomic distance 
is dictated by the nature of bonding. We expect that 
new phenomena will occur in these 'artificial molecules' 
(AMs) when the relative importance of Coulomb inter- 
action and single-particle tunneling is varied. Through 
the study of such phenomena, coupled dots may become 
a unique laboratory to explore electronic correlations. 

In this paper we analyze ground and excited few- 
particle states of realistic double quantum dots (DQDs) 
using a theoretical scheme that fully takes into account 
intra- and inter-dot many-body interactions. We focus 
on strongly confined vertically-coupled DQDs, and show 
that, for a given number of electrons, N , the ground 
state configuration undergoes non-trivial quantum tran- 
sitions as a function of the inter-dot distance d\ we iden- 
tify specific ranges of the DQD parameters characteriz- 
ing a 'coherent' molecular phase, where tunneling effects 



dominate, and a phase where the inter-dot interaction is 
purely electrostatic. At intermediate values of d we pre- 
dict a new phase that should be observable by means of 
transport experiments in the addition spectra and in the 
magnetization changes. 

From the theoretical point of view, a major difficulty 
is that we cannot make use of perturbative schemes 
PI in the calculation of the DQD many-body ground 
state, since we want to investigate all inter-dot coupling 
regimes. We thus write the many-body hamiltonian H 
within a generahzed Hubbard (GHH) approach [Q, and 
chose a basis set formed by suitable single-particle wave- 
functions localized on either dot {i = 1,2), and charac- 
terized by orbital quantum numbers, a (to be specified 
below), and by the spin quantum number, a. On this 
basis the many-body hamiltonian is 
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The first two addenda are the single-particle on-site and 
hopping term respectively: £„ are the single-particle en- 
ergies, t the tunneling parameter; h are the occupation 
numbers and c (ct) the creation (destruction) opera- 
tors. The third and fourth terms account for intra-dot 
Coulomb interaction between electrons with antiparallel 
and parallel spins, respectively: U and J are the intra- 
dot Coulomb and exchange integrals. Finally, the last 
term represents the inter-dot Coulomb coupling, U be- 
ing the inter-dot Coulomb integrals. These integrals are 
in turn expressed in terms of single-particle states: 
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with (/'ia(r) full three-dimensional single-particle wave- 
functions (obtained within the usual envelope-function 
approximation), e electronic charge and k^ relative di- 
electric constant. 



Here we shall consider a specific type of nanostruc- 
tures, namely vertically coupled cylindrical DQDs §; the 
ingredients entering the Hamiltonian can therefore be de- 
fined explicitly. For simplicity, V{r) can be assumed to 
be separable in the xy and z components; the profile is 
taken to be parabolic in the xy plane, with confinement 
energy Hlu, and a symmetric double quantum well (QW) 
along z ||5|,|6|. The eigenstates of the xy harmonic po- 
tential are the usual Fock-Darwin states |a) = \{n,m)) 
(n, m radial and angular quantum numbers) l^]. Along 
z, the QW thickness is such that only the lowest eigen- 
states (symmetric, |s), and antisymmetric \a)) are rele- 
vant for the low-energy spectrum. From these, we con- 
struct a complete set of states that are localized on ei- 
ther dot (see inset of Fig. 1): |1) = (|s) + \a)) /V2 and 
1 2) = (|s) — I a)) /V^. The basis we use is therefore the 
direct product \iaa) = \i) (g) |a) (g) |cr). The single-particle 
energies are Sa — £a + (^s + ea)/2, and the tunneling pa- 
rameter is i = (ea — es)/2, with e^, Eq double-well eigenen- 
ergies, and Sa oscillator energies. 

Using the above expressions, the three-dimensional 
Coulomb integrals and the tunneling parameter are com- 
puted directly from the single-particle states for each 
sample. Fig. 1 shows the result for t, Uaj3 and Uap 
[a^ j3= (0, 0)] calculated for a GaAs/AlAs DQD (k^ = 
12.98, effective mass 0.065 rrie) with fiuj = 10 meV, as a 
function of the interdot distance (barrier width) d. Note 
that already around d = 5 nm the inter-dot Coulomb 
integral exceeds the single-particle hopping parameter. 

To obtain the many-body energies and eigenstates, 
the Hamiltonian H is then diagonalized exactly for each 
value of N (total number of electrons) on each configu- 
ration subspace labeled by the quantum numbers S (z- 
component of the total spin) and M (z-component of the 
total angular momentum). 

This approach has two main advantages: First, it al- 
lows us to solve the many-body problem consistently in 
the different coupling regimes, from the limit where t 
dominates over the Coulomb integrals to the opposite 
limit where t is negligible. Second, it can provide quan- 
titative predictions for given DQD structures, since it 
contains no free parameters and uses realistic ingredients 
(i, Uaj3, Uaf3, Jafs) Calculated for each nanostructure 0. 

Fig. 2 shows the calculated ground-state energies E^f 
of correlated iV-particle states as a function of the inter- 
dot distance d [||. The three lowest excited states are 
also shown for comparison. As expected, when d is large, 
the system behaves as two isolated QDs. With decreas- 
ing d, some of the many-body excited states, favoured 
by Coulomb interactions, become lower in energy. For 
iV < 3 a single quantum phase transition occurs, below 
which the new ground state is a molecular-like state. For 
A^ > 3, two successive transitions take place, at d = db 
and d = da, and an intermediate non-trivial phase is pre- 
dicted to occur in the range da < d < db- Note that 
this phase is stable in a relatively large range of d values, 
which depends on the number of electrons. An accurate 



determination of da and db requires the correct inclusion 
of all inter-dot coupling terms, including the inter-dot 
Coulomb integrals. If the latter are neglected, all quan- 
tum transitions occur for smaller values of the inter-dot 
distance, and the d-range of the intermediate phase is un- 
derestimated significantly; had also the other many body 
terms been neglected, the intermediate phase would dis- 
appear leaving a simple crossover from a molecular to an 
atomic like regime at d = da — db. 

To understand the nature of the different phases it 
is useful to examine the many-body states in terms of 
the single Slater determinants that contribute to each 
of them. Here we discuss explicitly the 4-electron case 
with the help of the insets in Fig. 2, but the same rea- 
soning can be followed for the other cases. Both in the 
case of very small and very large interdot distances the 
ground state can be essentially described in terms of a 
single Slater determinant: for large values of d {d > db), 
the relevant configurations for the ground state have two 
electrons in the lowest level of each isolated dot; the \s) 
and I a) extended 'molecular' orbitals derived from the 
lowest 'atomic' states are of course almost degenerate 
and both filled with two electrons. In the opposite limit 
(d < da), by expanding the localized atomic orbitals in 
terms of molecular orbitals we recognize that the \s) state 
derived from the lowest atomic state is filled with two 
electrons, while the corresponding \a) molecular state is 
empty. The two remaining electrons occupy the next 
bonding molecular orbitals — derived from the higher p^ 
and py levels — with parallel spin, in such a way that 
S is maximized. This is the manifestation of Coulomb 
interaction which leads to Hund's rule for molecules. 

In these two extreme phases the single particle picture 
is essentially correct, provided that the appropriate basis 
set (either localized or extended orbitals) is used. In the 
intermediate phase, da < d < db, this is no longer true 
and the ground state is a mixture of different Slater de- 
terminants in any basis set H. In this sense again the 
intermediate phase exhibits an intrinsic many-body char- 
acter. Coulomb direct and exchange terms, responsible 
for the selection of the global quantum numbers, deter- 
mine a new ground state configuration, where both S and 
M are maximized. 

We obtain a clear evidence of the different electronic 
distribution in the three quantum phases by calculating 
the spin-dependent electronic pair correlation function, 
defined as 

5.^." {p, z', z") ^JdR (i>l (r')*^, (r")*." (r")*.' (r') 

where R = {p' + p")/2 and p = {p' - p"). Here 
p' = (a;', y'), p" = (x" ,y") are the in-plane spatial coordi- 
nates of the electrons in the pair, and z' and z" are their 
coordinates along z, that will be kept fixed at the center 
of either QD, i.e. in z\ or z^. r' = {p',z'), r" = {p",z"). 
a' , a" are the spin variables and can assume the values 
t or [. In Fig. 3 we plot for example g[,i{x, y, z' , z") for 
a double QD with A^ = 6 electrons, at three values of 



the inter-dot distance corresponding to different quan- 
tum phases. Here | represents the minority spins. For 
d < da (left column), it is indeed apparent that the pair 
correlation function is the same when both electrons are 
in the same dot [z' — z" = zi, panel (a)] or on differ- 
ent dots [z' — z\ and z" — Z2, panel (b)]. In this sense, 
the system behaves 'coherently'. For d > db (right col- 
umn), the maps indicate that the probability of finding 
two J,-electrons in the same dot is negligible, consistently 
with the picture of isolated dots. In the intermediate 
phase da < d < db this is no longer the case: the elec- 
tronic wavefunctions extend over both dots, and the g^^i 
pair correlation functions are very different depending on 
the location of both electrons in the same or in different 
dots. Note that transitions between different electronic 
configurations vs. d were recently identified theoretically 
also for classical coupled dots ||l^: We find that, for 
the small values of N considered here, the number of 
distinct phases and the spatial distribution of electrons 
(as reflected in their correlation functions) is drastically 
modified by quantum effects. 

The above findings are expected to be observable ex- 
perimentally. First, the calculated magnetic- field depen- 
dent addition spectra A{N) present clear signatures of 
the phase transitions described above, as illustrated in 
Fig. 4. Here A{N) = E{N) - E{N - 1) is obtained from 
the many-body ground state energies of the N- and (N- 
l)-electron systems, on the basis of single-particle states 
calculated in the presence of the external magnetic field 
B. The behaviour shown in the left, central, and right 
panels is representative of the three phases. Secondly, 
the changes in the magnetization induced by one elec- 
tron addition — also accessible experimentally |ll| — are 
expected to follow a different pattern in each phase (see 
the different sequence of quantum numbers in Fig. 2). 

A large experimental effort is presently devoted to 
transport experiments in double QDs. In most cases, 
the dots are obtained by gating a two-dimensional elec- 
tron gas (lateral confinement) , and their coupling can be 
tuned through a gate voltage [|l-15|. Indeed, a set of ex- 
periments has recently demonstrated a clear transition to 
a 'coherent' state with increasing coupling between dots 
|^,|5[. Outside this strong couphng regime, however, 
experiments performed in the lateral geometry have so- 
far evidenced classical interdot capacitance effects, prob- 
ably owing to the size of dots ||l6| . On the other hand, 
transport experiments are now available on DQD struc- 
tures with strong lateral confinement fabricated by com- 
bined growth and etching techniques ||5|,[l7|. The advan- 
tage is that the number of electrons in the structure is 
limited, while an accurate control on the inter-dot cou- 
pling is still possible by designing samples with appro- 
priate barrier thickness. Both aspects are important to 
enhance many-body effects and to explore intermediate 
coupling regimes. Recent experimental work on these 
structures has focused on the weak-coupling regimes f^ . 
We hope that further investigations will be stimulated by 
the present work, since the relevant transitions are now 



predicted quantitatively and should occur in a range of 
parameters that is accessible to state-of-the-art experi- 
ments. We expect that such studies will bring new in- 
sight into electron-electron interaction effects in coupled 
quantum nanostructures. 
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FIG. 1. Parameters entering the many-body Hamilto- 
nian of a double quantum dot with Tiuj = 10 meV, plotted 
as a function of the interdot distance d. The hopping coef- 
ficient t is shown together with two of the intra-dot and in- 
ter-dot Coulomb integrals. In the inset, the top panel shows 
the confinement potential V{z) (the barrier height is 200 
meV) and the corresponding symmetric and antisymmetric 
single-particle wavefunctions, \s) and \a); the bottom panel 
displays the localized states, |1) and |2), obtained as com- 
binations of \s) and \a), that are used as basis set for our 
calculation. 



FIG. 2. Energies of the ground state (thick line) and three 
lowest excited states (thin lines) as a function of the inter-dot 
distance, d, for a double QD occupied by A'^ electrons. For 
N = A, the prominent single particle configurations contribut- 
ing to the many-body ground state are shown in the insets. 
Note that in the intermediate phase a significant contribution 
comes also from other Slater determinants (see text). The -n 
molecular states are doubly degenerate because of the two- fold 
degeneracy of the second shell in the single-dots {px,Py)- 



FIG. 3. Electronic pair correlation functions for a double 
QD with N — 6 electrons, calculated at three values of the 
inter-dot distance, d, corresponding to the different quantum 
phases. The maps are plotted as a function of the difference 
between the in-plane coordinates of the electrons in the pair 
{x = x' — x" and y ~ y' — y"), while their z coordinates are 
fixed either at the center of the same dot [gi^i{x,y,Z\,z\); up- 
per panels] or at the center of different dots [gi,[{x, y, Zi, 22); 
lower panels]. Here J, represents the minority spins. 



FIG. 4. Cal- 

culated addition spectra, A{N) = E{N) - E{N - 1), as a 
function of magnetic field B for a double QD, calculated at 
three values of the inter-dot distance, d, corresponding to the 
different quantum phases. Labels indicate the ground state 
configuration of the A''-electron system. The energy zero is 
arbitrary. 




10 20 30 40 

inter-dot distance c/(nm) 



50 



M. Rontani et al.: Fig. 1 



761 

754 

747 

740 
611 

(D603 

^595 
■6)472 

(U466 

460 
342 

336 

330 




ground state 

three lowest excited states 



6 7 8 9 

inter-dot distance c/(nm) 



10 



M. Rontani et al.: Fig. 2 



d = 5 nm 



d = 7 nm 



d = 9 nm 



(a) 





(b) 






irr^ 



M. Rontani et al.: Fig. 3 



c/=0 nm ci=7 nm c/=9 nm 



30 




20 



> 

CD 

E 




^St lA- 



10 







'n,- ^A,+ 





- N=6 



N=5 



^.9 ^ 



N=4 



N=3 



0123012301234 

B(T) 



M. Rontani et al.: Fig. 4 



